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The promise of quantum computers is that certain computational tasks might be 
executed exponentially faster on a quantum processor than on a classical processor 1 . A 
fundamental challenge is to build a high-fidelity processor capable of running quantum 
algorithms in an exponentially large computational space. Here we report the use of a 
processor with programmable superconducting qubits 2 7 to create quantum states on 
53 qubits, corresponding to a computational state-space of dimension 2 53 (about 10 16 ). 
Measurements from repeated experiments sample the resulting probability 
distribution, which we verify using classical simulations. Our Sycamore processor takes 
about 200 seconds to sample one instance of a quantum circuit a million times-our 
benchmarks currently indicate that the equivalent task for a state-of-the-art classical 
supercomputer would take approximately 10,000 years. This dramatic increase in 
speed compared to all known classical algorithms is an experimental realization of 
quantum supremacy 814 for this specific computational task, heralding a much- 
anticipated computing paradigm. 


In the early 1980s, Richard Feynman proposed that a quantum computer 
would be an effective tool with which to solve problems in physics 
and chemistry, given that it is exponentially costly to simulate large 
quantum systems with classical computers 1 . Realizing Feynman's vision 
poses substantial experimental and theoretical challenges. First, can 
a quantum system be engineered to perform a computation in a large 
enough computational (Hilbert) space and with a low enough error 
rate to provide a quantum speedup? Second, can we formulate a prob¬ 
lem that is hard for a classical computer but easy for a quantum com¬ 
puter? By computing such a benchmark task on our superconducting 
qubit processor, we tackle both questions. Our experiment achieves 
quantum supremacy, a milestone on the path to full-scale quantum 
computing 814 . 


In reaching this milestone, we show that quantum speedup is achiev¬ 
able in a real-world system and is not precluded by any hidden physical 
laws. Quantum supremacy also heralds the era of noisy intermediate- 
scale quantum (NISQ) technologies 15 . The benchmark task we demon¬ 
strate has an immediate application in generating certifiable random 
numbers (S. Aaronson, manuscript in preparation); other initial uses 
for this new computational capability may include optimization 16,17 , 
machine learning 18-21 , materials science and chemistry 22-24 . However, 
realizing the full promise of quantum computing (using Shor's algorithm 
for factoring, for example) still requires technical leaps to engineer 
fault-tolerant logical qubits 25-29 . 

To achieve quantum supremacy, we made a number of techni¬ 
cal advances which also pave the way towards error correction. We 
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developed fast, high-fidelity gates that can be executed simultaneously 
across a two-dimensional qubit array. We calibrated and benchmarked 
the processor at both the component and system level using a powerful 
new tool: cross-entropy benchmarking 11 . Finally, we used component- 
level fidelities to accurately predict the performance of the whole sys¬ 
tem, further showing that quantum information behaves as expected 
when scaling to large systems. 

A suitable computational task 

To demonstrate quantum supremacy, we compare our quantum proces¬ 
sor against state-of-the-art classical computers in the task of sampling 
the output of a pseudo-random quantum circuit 111314 . Random circuits 
are a suitable choice for benchmarking because they do not possess 
structure and therefore allow for limited guarantees of computational 
hardness 1012 . We design the circuits to entangle a set of quantum bits 
(qubits) by repeated application of single-qubit and two-qubit logi¬ 
cal operations. Sampling the quantum circuit's output produces a set 
of bitstrings, for example {0000101,1011100,...}. Owing to quantum 
interference, the probability distribution of the bitstrings resembles 
a speckled intensity pattern produced by light interference in laser 
scatter, such that some bitstrings are much more likely to occur than 
others. Classically computing this probability distribution becomes 
exponentially more difficult as the number of qubits (width) and number 
of gate cycles (depth) grow. 

We verify that the quantum processor is working properly using a 
method called cross-entropy benchmarking 111214 , which compares how 
often each bitstring is observed experimentally with its corresponding 
ideal probability computed via simulation on a classical computer. For 
a given circuit, we collect the measured bitstrings {xj and compute the 
linear cross-entropy benchmarking fidelity 111314 (see also Supplementary 
Information), which is the mean of the simulated probabilities of the 
bitstrings we measured: 


•^XEB -2 n (P(X i )) i 1 (1) 

where n is the number of qubits, P(x t ) is the probability of bitstring x z 
computed for the ideal quantum circuit, and the average is over the 
observed bitstrings. Intuitively, :r XEB is correlated with how often we 
sample high-probability bitstrings. When there are no errors in the 
quantum circuit, the distribution of probabilities is exponential (see 
Supplementary Information), and sampling from this distribution will 
produce :F xeb = 1. On the other hand, sampling from the uniform 
distribution will give (P(x))i = 1/2" and produce J^ XEB = 0. Values of J^ XEB 
between 0 and 1 correspond to the probability that no error has occurred 
while running the circuit. The probabilities P(x t ) must be obtained from 
classically simulating the quantum circuit, and thus computing J^ XEB is 
intractable in the regime of quantum supremacy. Flowever, with certain 
circuit simplifications, we can obtain quantitative fidelity estimates of 
a fully operating processor running wide and deep quantum circuits. 

Our goal is to achieve a high enough J^ XEB for a circuit with sufficient 
width and depth such that the classical computing cost is prohibitively 
large. This is a difficult task because our logic gates are imperfect and 
the quantum states we intend to create are sensitive to errors. A single 
bit or phase flip over the course of the algorithm will completely shuffle 
the speckle pattern and result in close to zero fidelity 11 (see also Sup¬ 
plementary Information). Therefore, in order to claim quantum suprem¬ 
acy we need a quantum processor that executes the program with 
sufficiently low error rates. 

Building a high-fidelity processor 

We designed a quantum processor named 'Sycamore' which consists 
of a two-dimensional array of 54 transmon qubits, where each qubit is 
tunably coupled to four nearest neighbours, in a rectangular lattice. The 
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Fig. 1 1 The Sycamore processor, a. Layout of processor, showing a rectangular 
array of 54 qubits (grey), each connected to its four nearest neighbours with 
couplers (blue). The inoperable qubit is outlined, b. Photograph of the 
Sycamore chip. 


connectivity was chosen to be forward-compatible with error correc¬ 
tion using the surface code 26 . A key systems engineering advance of this 
device is achieving high-fidelity single- and two-qubit operations, not 
just in isolation but also while performing a realistic computation with 
simultaneous gate operations on many qubits. We discuss the highlights 
below; see also the Supplementary Information. 

In a superconducting circuit, conduction electrons condense into a 
macroscopic quantum state, such that currents and voltages behave 
quantum mechanically 2,30 . Our processor uses transmon qubits 6 , which 
can be thought of as nonlinear superconducting resonators at 5-7 GFIz. 
The qubit is encoded as the two lowest quantum eigenstates of the 
resonant circuit. Each transmon has two controls: a microwave drive 
to excite the qubit, and a magnetic flux control to tune the frequency. 
Each qubit is connected to a linear resonator used to read out the qubit 
state 5 . As shown in Fig. 1, each qubit is also connected to its neighbouring 
qubits using a new adjustable coupler 31,32 . Our coupler design allows us 
to quickly tune the qubit-qubit coupling from completely off to 40 MFIz. 
One qubit did not function properly, so the device uses 53 qubits and 
86 couplers. 

The processor is fabricated using aluminium for metallization and 
Josephson junctions, and indium for bump-bonds between two silicon 
wafers. The chip is wire-bonded to a superconducting circuit board 
and cooled to below 20 mK in a dilution refrigerator to reduce ambient 
thermal energy to well below the qubit energy. The processor is con¬ 
nected through filters and attenuators to room-temperature electronics. 



















which synthesize the control signals. The state of all qubits can be read 
simultaneously by using a frequency-multiplexing technique 33 34 . We use 
two stages of cryogenic amplifiers to boost the signal, which is digitized 
(8 bits at 1 GHz) and demultiplexed digitally at room temperature. In 
total, we orchestrate 277 digital-to-analog converters (14 bits at 1 GHz) 
for complete control of the quantum processor. 

We execute single-qubit gates by driving 25-ns microwave pulses reso¬ 
nant with the qubit frequency while the qubit-qubit coupling is turned 
off. The pulses are shaped to minimize transitions to higher transmon 
states 35 . Gate performance varies strongly with frequency owing to two- 
level-system defects 36 37 , stray microwave modes, coupling to control 
lines and the readout resonator, residual stray coupling between qubits, 
flux noise and pulse distortions. We therefore optimize the single-qubit 
operation frequencies to mitigate these error mechanisms. 

We benchmark single-qubit gate performance by using the cross¬ 
entropy benchmarking protocol described above, reduced to the single¬ 
qubit level (n = 1), to measure the probability of an error occurring 
during a single-qubit gate. On each qubit, we apply a variable number 
m of randomly selected gates and measure J^ XEB averaged over many 
sequences; as m increases, errors accumulate and average J^ XEB decays. 
We model this decay by [1 - ej(\ - 1 /D 2 )] m where e x is the Pauli error prob¬ 
ability. The state (Hilbert) space dimension term, D = 2 n , which equals 
2 for this case, corrects for the depolarizing model where states with 
errors partially overlap with the ideal state. This procedure is similar to 
the more typical technique of randomized benchmarking 27 ' 38 ' 39 , but 
supports non-CIifford-gate sets 40 and can separate out decoherence 
error from coherent control error. We then repeat the experiment with 
all qubits executing single-qubit gates simultaneously (Fig. 2), which 
shows only a small increase in the error probabilities, demonstrating 
that our device has low microwave crosstalk. 

We perform two-qubit iSWAP-like entangling gates by bringing neigh¬ 
bouring qubits on-resonance and turning on a 20-MHz coupling for 12 ns, 
which allows the qubits to swap excitations. During this time, the qubits 
also experience a controlled-phase (CZ) interaction, which originates 
from the higher levels of the transmon. The two-qubit gate frequency 
trajectories of each pair of qubits are optimized to mitigate the same error 
mechanisms considered in optimizing single-qubit operation frequencies. 

To characterize and benchmarkthe two-qubit gates, we run two-qubit 
circuits with m cycles, where each cycle contains a randomly chosen 
single-qubit gate on each of the two qubits followed by a fixed two-qubit 
gate. We learn the parameters of the two-qubit unitary (such as the 
amount of iSWAP and CZ interaction) by using J^ XEB as a cost function. 
After this optimization, we extract the per-cycle error e lc from the decay 
of J^xeb with m, and isolate the two-qubit error e 2 by subtracting the two 
single-qubit errors e x . We find an average e 2 of 0.36%. Additionally, we 
repeat the same procedure while simultaneously running two-qubit 
circuits for the entire array. After updating the unitary parameters to 
account for effects such as dispersive shifts and crosstalk, we find an 
average e 2 of 0.62%. 

For the full experiment, we generate quantum circuits using the two- 
qubit unitaries measured for each pair during simultaneous operation, 
rather than a standard gate for all pairs. The typical two-qubit gate is a 
full iSWAP with l/6th of a full CZ. Using individually calibrated gates in 
no way limits the universality of the demonstration. One can compose, 
for example, controlled-NOT (CNOT) gates from 1-qubit gates and two 
of the unique 2-qubit gates of any given pair. The implementation of 
high-fidelity Textbook gates 7 natively, such as CZ or ^iSWAP, is work 
in progress. 

Finally, we benchmark qubit readout using standard dispersive meas¬ 
urement 41 . Measurement errors averaged over the 0 and 1 states are 
shown in Fig. 2a. We have also measured the error when operating all 
qubits simultaneously, by randomly preparing each qubit in the 0 or 1 
state and then measuring all qubits for the probability of the correct 
result. We find that simultaneous readout incurs only a modest increase 
in per-qubit measurement errors. 
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Fig. 21 System-wide Pauli and measurement errors, a. Integrated histogram 
(empirical cumulative distribution function, ECDF) of Pauli errors (black, green, 
blue) and readout errors (orange), measured on qubits in isolation (dotted lines) 
and when operating all qubits simultaneously (solid). The median of each 
distribution occurs at 0.50 on the vertical axis. Average (mean) values are shown 
below, b, Heat map showing single- and two-qubit Pauli errors e x (crosses) and e 2 
(bars) positioned in the layout of the processor. Values are shown for all qubits 
operating simultaneously. 


Having found the error rates of the individual gates and readout, we 
can model the fidelity of a quantum circuit as the product of the prob¬ 
abilities of error-free operation of all gates and measurements. Our 
largest random quantum circuits have 53 qubits, 1,113 single-qubit gates, 
430 two-qubit gates, and a measurement on each qubit, for which we 
predict a total fidelity of 0.2%. This fidelity should be resolvable with a 
few million measurements, since the uncertainty on ^ XEB is V^N S > where 
N s is the number of samples. Our model assumes that entangling larger 
and larger systems does not introduce additional error sources beyond 
the errors we measure at the single- and two-qubit level. In the next 
section we will see how well this hypothesis holds up. 

Fidelity estimation in the supremacy regime 

The gate sequence for our pseudo-random quantum circuit generation 
is shown in Fig. 3. One cycle of the algorithm consists of applying 
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Cycle 1 


Fig. 31 Control operations for the quantum supremacy circuits, a. Example 
quantum circuit instance used in our experiment. Every cycle includes a layer 
each of single- and two-qubit gates. The single-qubit gates are chosen randomly 
from {4x, 4Y, y[W], where W=(X+ Y)/42 and gates do not repeat sequentially. 
The sequence of two-qubit gates is chosen according to a tiling pattern, 
coupling each qubit sequentially to its four nearest-neighbour qubits. The 


couplers are divided into four subsets (ABCD), each of which is executed 
simultaneously across the entire array corresponding to shaded colours. Here 
we show an intractable sequence (repeat ABCDCDAB); we also use different 
coupler subsets along with a simplifiable sequence (repeat EFGHEFGH, not 
shown) that can be simulated on a classical computer, b. Waveform of control 
signals for single- and two-qubit gates. 


single-qubit gates chosen randomly from {4x, 4Y, 4W} on all qubits, 
followed by two-qubit gates on pairs of qubits. The sequences of gates 
which form the 'supremacy circuits' are designed to minimize the circuit 
depth required to create a highly entangled state, which is needed for 
computational complexity and classical hardness. 

Although we cannot compute :F xeb in the supremacy regime, we can 
estimate it using three variations to reduce the complexity of the circuits. 
In 'patch circuits', we remove a slice of two-qubit gates (a small fraction 
of the total number of two-qubit gates), splitting the circuit into two 
spatially isolated, non-interacting patches of qubits. We then compute 
the total fidelity as the product of the patch fidelities, each of which can 
be easily calculated. In 'elided circuits', we remove only a fraction of the 
initial two-qubit gates along the slice, allowing for entanglement 
between patches, which more closely mimics the full experiment while 
still maintaining simulation feasibility. Finally, we can also run full 
'verification circuits', with the same gate counts as our supremacy cir¬ 
cuits, but with a different pattern for the sequence of two-qubit gates, 
which is much easier to simulate classically (see also Supplementary 
Information). Comparison between these three variations allows us to 
track the system fidelity as we approach the supremacy regime. 

We first check that the patch and elided versions of the verification 
circuits produce the same fidelity as the full verification circuits up to 
53 qubits, as shown in Fig. 4a. For each data point, we typically collect 
A/ s = 5 x 10 6 total samples over ten circuit instances, where instances 
differ only in the choices of single-qubit gates in each cycle. We also 
show predicted J^ XEB values, computed by multiplying the no-error prob¬ 
abilities of single- and two-qubit gates and measurement (see also Sup¬ 
plementary Information). The predicted, patch and elided fidelities all 
show good agreement with the fidelities of the corresponding full cir¬ 
cuits, despite the vast differences in computational complexity and 
entanglement. This gives us confidence that elided circuits can be used 
to accurately estimate the fidelity of more-complex circuits. 

The largest circuits for which the fidelity can still be directly verified 
have 53 qubits and a simplified gate arrangement. Performing random 
circuit sampling on these at 0.8% fidelity takes one million cores 130 
seconds, corresponding to a million-fold speedup of the quantum pro¬ 
cessor relative to a single core. 

We proceed now to benchmark our computationally most difficult 
circuits, which are simply a rearrangement of the two-qubit gates. In 
Fig. 4b, we show the measured :F xeb for 53-qubit patch and elided ver¬ 
sions of the full supremacy circuits with increasing depth. For the larg¬ 
est circuit with 53 qubits and 20 cycles, we collected N s = 30 x 10 6 samples 
over ten circuit instances, obtaining jr XEB = (2.24 ±0.21) x 10 -3 for the 
elided circuits. With 5a confidence, we assert that the average fidelity 


of running these circuits on the quantum processor is greater than at 
least 0.1%. We expect that the full data for Fig. 4b should have similar 
fidelities, but since the simulation times (red numbers) take too long to 
check, we have archived the data (see 'Data availability' section). The 
data is thus in the quantum supremacy regime. 

The classical computational cost 

We simulate the quantum circuits used in the experiment on classical 
computers for two purposes: (1) verifying our quantum processor and 
benchmarking methods by computing J^ XEB where possible using sim¬ 
plifiable circuits (Fig. 4a), and (2) estimating J^ XEB as well as the classical 
cost of sampling our hardest circuits (Fig. 4b). Up to 43 qubits, we use 
a Schrodinger algorithm, which simulates the evolution of the full quan¬ 
tum state; thejulich supercomputer (with 100,000 cores, 250 terabytes) 
runs the largest cases. Above this size, there is not enough random access 
memory (RAM) to store the quantum state 42 . For larger qubit numbers, 
we use a hybrid Schrodinger-Feynman algorithm 43 running on Google 
data centres to compute the amplitudes of individual bitstrings. This 
algorithm breaks the circuit up into two patches of qubits and efficiently 
simulates each patch using a Schrodinger method, before connecting 
them using an approach reminiscent of the Feynman path-integral. 
Although it is more memory-efficient, the Schrodinger-Feynman algo¬ 
rithm becomes exponentially more computationally expensive with 
increasing circuit depth owing to the exponential growth of paths with 
the number of gates connecting the patches. 

To estimate the classical computational cost of the supremacy circuits 
(grey numbers in Fig. 4b), we ran portions of the quantum circuit simu¬ 
lation on both the Summit supercomputer as well as on Google clusters 
and extrapolated to the full cost. In this extrapolation, we account for 
the computation cost of sampling by scaling the verification cost with 
JF XEB , for example 43 ' 44 , a 0.1% fidelity decreases the cost by about 1,000. 
On the Summit supercomputer, which is currently the most powerful 
in the world, we used a method inspired by Feynman path-integrals that 
is most efficient at low depth 44-47 . At m = 20 the tensors do not reason¬ 
ably fit into node memory, so we can only measure runtimes up to m= 14, 
for which we estimate that sampling three million bitstrings with 1% 
fidelity would require a year. 

On Google Cloud servers, we estimate that performing the same task 
for m =20 with 0.1% fidelity using the Schrodinger-Feynman algorithm 
would cost 50 trillion core-hours and consume one petawatt hour of 
energy. To put this in perspective, it took 600 seconds to sample the 
circuit on the quantum processor three million times, where sampling 
time is limited by control hardware communications; in fact, the net 
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Fig. 41 Demonstrating quantum supremacy, a. Verification of benchmarking 
methods. ^ XEB values for patch, elided and full verification circuits are 
calculated from measured bitstrings and the corresponding probabilities 
predicted by classical simulation. Here, the two-qubit gates are applied in a 
simplifiable tiling and sequence such that the full circuits can be simulated out 
to n = 53, m = 14 in a reasonable amount of time. Each data point is an average over 
ten distinct quantum circuit instances that differ in their single-qubit gates (for n 
= 39,42 and 43 only two instances were simulated). For each n, each instance is 
sampled with /V s of 0.5-2.5 million. The black line shows the predicted j^ XEB based 
on single- and two-qubit gate and measurement errors. The close 
correspondence between all four curves, despite their vast differences in 


complexity, justifies the use of elided circuits to estimate fidelity in the 
supremacy regime, b. Estimating ^ XEB in the quantum supremacy regime. Here, 
the two-qubit gates are applied in a non-simplifiable tiling and sequence for 
which it is much harder to simulate. For the largest elided data (n = 53, m = 20, 
total N s = 30 million), we find an average jf xeb > 0.1% with 5aconfidence, where o 
includes both systematic and statistical uncertainties. The corresponding full 
circuit data, not simulated but archived, is expected to show similarly 
statistically significant fidelity. For m = 20, obtaining a million samples on the 
quantum processortakes200 seconds, whereasan equal-fidelity classical 
sampling would take 10,000 years on a million cores, and verifying the fidelity 
would take millions of years. 


quantum processor time is only about 30 seconds. The bitstring samples 
from all circuits have been archived online (see 'Data availability' section) 
to encourage development and testing of more advanced verification 
algorithms. 

One may wonder to what extent algorithmic innovation can enhance 
classical simulations. Our assumption, based on insights from complex¬ 
ity theory 1113 , is that the cost of this algorithmic task is exponential in 
circuit size. Indeed, simulation methods have improved steadily over the 
past few years 42 50 . We expect that lower simulation costs than reported 
here will eventually be achieved, but we also expect that they will be 
consistently outpaced by hardware improvements on larger quantum 
processors. 

Verifying the digital error model 

A key assumption underlying the theory of quantum error correction 
is that quantum state errors may be considered digitized and local¬ 
ized 38,51 . Under such a digital model, all errors in the evolving quantum 
state may be characterized by a set of localized Pauli errors (bit-flips or 
phase-flips) interspersed into the circuit. Since continuous amplitudes 
are fundamental to quantum mechanics, it needs to be tested whether 
errors in a quantum system could be treated as discrete and probabil¬ 
istic. Indeed, our experimental observations support the validity of 
this model for our processor. Our system fidelity is well predicted by a 
simple model in which the individually characterized fidelities of each 
gate are multiplied together (Fig. 4). 

To be successfully described by a digitized error model, a system 
should be low in correlated errors. We achieve this in our experiment by 


choosing circuits that randomize and decorrelate errors, by optimizing 
control to minimize systematic errors and leakage, and by designing 
gates that operate much faster than correlated noise sources, such as 
l//flux noise 37 . Demonstrating a predictive uncorrelated error model 
up to a Hilbert space of size 2 53 shows that we can build a system where 
quantum resources, such as entanglement, are not prohibitively fragile. 

The future 

Quantum processors based on superconducting qubits can now perform 
computations in a Hilbert space of dimension 2 53 » 9 x io 15 , beyond the 
reach of the fastest classical supercomputers available today. To our 
knowledge, this experiment marks the first computation that can be 
performed only on a quantum processor. Quantum processors have 
thus reached the regime of quantum supremacy. We expect that their 
computational power will continue to grow at a double-exponential 
rate: the classical cost of simulating a quantum circuit increases expo¬ 
nentially with computational volume, and hardware improvements will 
probably follow a quantum-processor equivalent of Moore's law 52,53 , 
doubling this computational volume every few years. To sustain the 
double-exponential growth rate and to eventually offer the computa¬ 
tional volume needed to run well known quantum algorithms, such as 
the Shor or Grover algorithms 25,54 , the engineering of quantum error 
correction will need to become a focus of attention. 

The extended Church-Turing thesis formulated by Bernstein and 
Vazirani 55 asserts that any 'reasonable' model of computation can be 
efficiently simulated by a Turing machine. Our experiment suggests 
that a model of computation may now be available that violates this 
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assertion. We have performed random quantum circuit sampling in 
polynomial time using a physically realizable quantum processor (with 
sufficiently low error rates), yet no efficient method is known to exist for 
classical computing machinery. As a result of these developments, quan¬ 
tum computing is transitioning from a research topic to a technology 
that unlocks new computational capabilities. We are only one creative 
algorithm away from valuable near-term applications. 

Data availability 

The datasets generated and analysed for this study are available at our 
public Dryad repository (https://doi.org/10.5061/dryad.k6tlrj8). 
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